Work functions of self-assembled monolayers on metal surfaces 



Paul C. Rush and Geert Brocks 
Computational Materials Science, Faculty of Science and Technology and MESA+ Institute for Nanotechnology, 
University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands 
(Dated: February 6, 2008) 

Using first-principles calculations we show that the work function of noble metals can be decreased 
or increased by up to 2 eV upon the adsorption of self-assembled monolayers of organic molecules. 
We identify the contributions to these changes for several (fluorinated) thiolate molecules adsorbed 
on Ag(lll), Au(lll) and Pt(lll) surfaces. The work function of the clean metal surfaces increases 
in this order, but adsorption of the monolayers reverses the order completely. Bonds between the 
thiolate molecules and the metal surfaces generate an interface dipole, whose size is a function of 
the metal, but it is relatively independent of the molecules. The molecular and bond dipoles can 
then be added to determine the overall work function. 

PACS numbers: 73.30.-(-y, 73.61.-Ph, 68.43.-h 



Recent advances in molecular electronics, where or- 
ganic molecules constitute active materials in electronic 
devices, have created a large interest in metal organic in- 
terfaces Q . Transport of charge carriers across the inter- 
faces between metal electrodes and the organic material 
often determines the performance of a device 0. Or- 
ganic semiconductors differ from inorganic ones as they 
are composed of molecules and intermolecular forces are 
relatively weak. In a bulk material this increases the 
importance of electron-phonon and electron-electron in- 
teractions 13 . At a metal organic interface the energy 
barrier for charge carrier injection into the organic mate- 
rial is often determined by the formation of an interface 
dipole localized at the first molecular layer. The inter- 
face dipole can be extracted by monitoring the change 
in the metal surface work function after deposition of an 
organic layer 1, 4J. 

Atoms and molecules that are physisorbed on a metal 
surface usually decrease the work function, as the Pauli 
repulsion between the molecular and surface electrons 
decreases the surface dipole 0, 0|. Chemisorption can 
give an increase or a decrease of the work function, and 
can even lead to counterintuitive results 0, Self- 
assembled monolayers (SAMs) are exemplary systems to 
study the effect of chemisorbed organic molecules upon 
metal work functions . More specifically, alkyl thiolate 
(C„H2„+iS) SAMs on the gold (111) surface are ainong 
the most extensively studied systems [l^ [ill, [l^, llM ll3 ■ 
The sulphur atoms of the thiolate molecules form stable 
bonds to the gold surface and their alkyl tails are close 
packed, which results in a well ordered monolayer. SAMs 
with similar structures are formed by alkyl thiolates on 
a range of other (noble) metal surfaces 10, 14, 15]. 

Often the change in work function upon adsorption 
of a SAM is interpreted mainly in terms of the dipole 
moments of the individual thiolate molecules, whereas 
only a minor role is attributed to the change induced 
by chemisorption [ll| This assumption turns 

out to be reasonable for adsorption of methyl thiolate 



(CHaS) on Au(lll) fl3(, but for CH3S on Cu(lll) it is 
not |l4| . In this paper we apply first- principles calcula- 
tions to study the interface dipoles and the work function 
change induced by adsorption of thiolate SAMs. 

In particular, we analyze the contributions of 
chemisorption and of the molecular dipoles to uncover 
the effects of charge reordering at the interface. The 
chemical bonds between the thiolate molecules and the 
metal surfaces generate an interface dipole. We find 
that this dipole strongly depends upon the metal, but 
it is nearly independent of the electronegativity of the 
molecules. The size and direction of the interface dipole 
are such that it overcompensates for the difference be- 
tween the clean metal work functions. This results in the 
SAM adsorbed on the highest work function metal hav- 
ing the lowest work function and vice versa. Modifying 
the molecular tails allows one to vary the absolute size of 
the work function over a range of more than 2 eV. 

Since alkyl thiolate molecules form SAMs with a sim- 
ilar structure on (111) surfaces of several noble metals, 
they are ideal model systems for studying metal organic 
interfaces. By varying the relative electronegativity of 
surface and molecules one can induce electron transfer 
and create an interface dipole, without completely re- 
arranging the interface structure. The electronegativity 
of a metal substrate is given by its work function. We 
consider the (111) surfaces of three metals that have a 
substantially different work function, but the same crys- 
tal structure and a similar lattice parameter: Ag (4.5 eV, 
2.89A ), Au (5.3 eV, 2.88A ) and Pt (6.1 eV, 2.77A ). 

One would also like to vary the molecule's electroneg- 
ativity without changing the structure of the SAM. This 
can be achieved by fluorinating the alkyl tails of thiolate 
molecules, which increases their electronegativity (lOi] . 
However, fluorinating the alkyl tails also reverses the po- 
larity of the thiolate molecules and one has to separate 
this electrostatic effect from the charge reordering caused 
by chemisorption. In this paper we study the short chain 
thiolates CH3S, C2H5S, CF3S, and CF3CH2S. 
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Density functional theory (DFT) calculations are car- 
ried out using the projector augmented wave (PAW) 
method 0, 0|, a plane wave basis set and the PW91 
generalized gradient approximation (GGA) functional, as 
implemented in the VASP program 0, ^ . We use su- 
percells containing a slab of at least five layers of metal 
atoms with a SAM adsorbed on one side of the slab and 
a vacuum region of ~ 12 A. The Brillouin zone of the 
{VS X \/3)R.30° surface unit cell is sampled by a 11 x 11 
k-point grid. The plane wave kinetic energy cutoff is 450 
eV. To avoid interactions between periodic images of the 
slab we apply a dipole correction |20j. The geometry of 
the SAM is optimized, as well as the positions of the top 
two layers of metal atoms. The atoms in the remaining 
metal layers are fixed at their bulk positions. The opti- 
mized bulk lattice parameters are 2.93, 2.94 and 2.79 A 
for Ag, Au and Pt, respectively. 

The work function is given hy W — V{oo) — Ep, where 
V{oo) is the asymptotic electrostatic potential in vac- 
uum, and Ep is the Fermi energy of the bulk metal. 
V{oo) is extracted from the plane averaged potential 
V{z) = JJ^V{x,y,z)dxdy, with A the area of the 
surface unit cell. In practice, V{z) reaches an asymp- 
totic value within a distance of 5 A from the surface. 
Accurate values of the Fermi energy are obtained follow- 
ing the procedure outlined in Ref. I21I By varying the 
computational parameters discussed above we estimate 
that the work functions are converged to within 0.05 eV. 
Typically DFT calculations give work functions that are 
within ^ 0.1-0.2 cV of the experimental values, although 
occasionally somewhat larger deviations are found. 

The (V3 X V3)R30° structure of CH3S on Au(lll) 
has been studied in several first-principles calculations 
[IIIIII2II2II23. We find basically the same optimized 
geometry as obtained in those calculations. Several struc- 
tures exist that have a slightly different geometry, but are 
very close in energy, such as a c(4 x 2) superstructure . 
We find that the work functions of these structures are 
within 0.1 eV of that of the simpler structure, so we will 
not discuss these superstructures here. 

The (v^ X •\/3)R30° structure is also a good starting 
point for studying other systems. Thiolates with longer 
alkyl tails on Au(lll) adopt this structure, as does CH3S 
on Pt(lll), as well as alkyl thiolates on Au(lll) whose 
end groups are fluorinated 0, 0| . Thiolates with long 
alkyl tails on Ag(lll) form a somewhat denser packing, 
whereas long fluorinated alkyl thiolates form a somewhat 
less dense packing 10]. To analyze the work function we 
use optimized (\/3x-\/3)R30° structures for all our SAMs. 
We find that varying the packing density only introduces 
a scaling factor to the work function change 0| . 

Tabic m lists the calculated work functions. The work 
functions of the clean Au and Ag surfaces agree with the 
experimental values [l^ll^, but that of Pt is ~ 0.3 eV 
too low [23. The latter can be attributed to the GGA 
functional. Using the local density approximation (LDA) 



TABLE I: Calculated work functions W (eV) of clean (111) 
surfaces and of surfaces covered by SAMs in a (\/3x^)R30° 
structure. 





clean 


CH3S 


C2H5S 


CF3S 


CF3CH2S 


Ag 


4.50 


3.95 


4.13 


6.14 


6.30 


Au 


5.25 


3.81 


3.93 


5.97 


6.27 


Pt 


5.84(6.14") 


3.45 


3.47 


5.68 


5.87 



'LDA value 



the calculated work function of Pt(lll) is 6.14 eV, which 
agrees with experiment. In other cases the difference be- 
tween the work functions calculated with GGA and LDA 
functionals is much smaller. For instance, the GGA and 
LDA work functions of the SAMs on Pt are within 0.02 
eV of one another. We will use the GGA values through- 
out this paper. The trend in the work functions of the 
SAM covered surfaces agrees well with experimental ob- 
servations j2[, 11.. JL2] . The experimental work function 
shifts with respect to the clean surface are sometimes 
somewhat smaller than the calculated ones [isj . 

The first observation one can make by comparing the 
numbers in Table within columns is that on SAM cov- 
ered surfaces the work function decreases in the order 
Ag, Au, Pt. This is striking, since the work function 
of the clean metal surfaces clearly increases in this or- 
der. Secondly, comparing the numbers within rows one 
finds that the work functions of the fluorinated alkyl thi- 
olate covered surfaces are 2-2.5 eV higher than of the 
non- fluorinated ones. We will argue that the first ob- 
servation can be ascribed to the interface dipole formed 
upon chemisorption. This interface dipole is independent 
of the molecular tails. The second observation will be in- 
terpreted in terms of the individual molecular dipoles. 

In order to visualize the charge reordering at the sur- 
face upon adsorption of the SAM, we calculate the differ- 
ence electron density An. It is obtained by subtracting 
from the total electron density rttot of the SAM on the 
surface, the electron density risurf of the clean surface and 
that of the free standing SAM tt-sam- '^smf and nsAM are 
obtained in two separate calculations of a clean surface 
and a free standing SAM, respectively, with their struc- 
tures frozen in the adsorbed geometry. As an example, 
Fig.CJshows An for SAMs of CF3S and CH3S on Ag(lll). 

Fig. n illustrates that An is localized mainly at the 
metal-SAM interface, i.e. near the sulphur atoms and 
the metal atoms in the first surface layers. In case of ad- 
sorption on Ag, electrons are transferred from the metal 
to the molecule, which results in an increase of the elec- 
tron density on the sulphur atoms and a decrease on the 
surface metal atoms. The charge transfer does not de- 
pend strongly on the molecule, compare Figs. n[^a,b) to 
(c,d). This is somewhat surprising, since the electroneg- 
ativity of CF3S is much higher than that of CH3S. 

Very often a charge transfer between two systems is in- 
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TABLE II: Dipole per molecule A/i, from the change in work function upon adsorption. The (perpendicular) molecular dipole 
moment ^sam in a free standing SAM. The chemisorption dipole moment is ^chcm = A/i — fisAM- All values are in D. 





CHsS 


C2H5S 


Ag 
CF3S 


CF3CH2S 


CH3S 


C2H5S 


Au 
CF3S 


CF3CH2S 


CH3S 


C2H5S 


Pt 
CF3S 


CF3CH2S 




-0.32 


-0.22 


0.97 


1.07 


-0.86 


-0.79 


0.43 


0.61 


-1.28 


-1.27 


-0.08 


0.02 


MSAM 


-0.88 


-0.79 


0.44 


0.50 


-0.88 


-0.81 


0.44 


0.53 


-0.86 


-0.80 


0.37 


0.47 


/^chcm 


0.56 


0.57 


0.53 


0.57 


0.02 


0.02 


-0.01 


0.08 


-0.42 


-0.47 


-0.45 


-0.45 




-0.004 0.000 0.004 



FIG. 1: (Color online) Difference electron density An — 
n-tot — risuri — nsAM for CF3S on Ag(lll), (a) as function of 
z, averaged over the xy plane, in units of A ^; (b) as an 
isodensity surface; (c), (d) the same for CH3S on Ag(lll). 



terpreted in terms of their relative electronegativity. For 
a metal surface the latter is simply the work function 
Wciean- For a molecule the Mulliken electronegativity 
Xm is defined as the average of the ionization poten- 
tial and the electron affinity and considered to be the 
molecular equivalent of a chemical potential |2^. We 
find Xm = 5.4 eV for the CH3S and CH3CH2S molecules. 
Since xm is close to Wdean for Au(lll), this would ex- 
plain the lack of electron transfer upon adsorption of 
these molecules However, the calculated xm for 

CF3S and CF3CH2S are much higher, i.e. 6.9 eV and 6.1 
eV, respectively. Yet this does not result in a markedly 
increased electron transfer to these molecules, as Fig. Q] 
indicates. It means that xm is not a generally suitable 
parameter to predict the amount of charge transfer be- 



tween surface and molecules, xm reflects the relative 
stability of charged molecular states. In particular, for 
the thiolates xm reflects the ability of the (fluorinated) 
alkyl chains to stabilize or screen charge that resides on 
the sulphur atom. We suggest that this is not important 
in case of adsorbed molecules, as the metal surface takes 
over this role. 

Meanwhile, Fig. ^ suggests the following analysis. 
From the change in the work function upon adsorption 
of the SAM, AW = W - Wcican, see Table |l| one can 
obtain the change of the surface dipole upon adsorption, 
A/x = eoAAW/e (with eg the permittivity of vacuum and 
A the area of the surface unit cell). Since the unit cell 
contains one molecule, A/i is the change in the surface 
dipole per adsorbed molecule. The results are shown in 
Table Hll A/i contains contributions from the charge re- 
ordering at the interface due to chemisorption, as well as 
from the dipole moments of the individual molecules. 

The latter can be accounted for by calculating the 
dipole moment /^sam per molecule of free standing SAMs, 
i.e. without the presence of a metal surface. We fo- 
cus upon the component of the dipole that is perpen- 
dicular to the surface, since the other components do 
not contribute to the work function. As the calculation 
uses a full monolayer of molecules, it incorporates the ef- 
fect on each molecule of the depolarizing fleld caused by 
the dipoles of all surrounding molecules. The calculated 
A*SAM are given in Table Ull The structure of a SAM is 
flxed in its adsorption geometry, which is similar for the 
three metal surfaces. Therefore, the fJ-sAM values for ad- 
sorption on Ag, Au, and Pt in Tablelllldiffer only slightly. 
Of course /^sam depends upon the molecule. In CH3S 
and CH3CH2S the dipole points from the sulphur atom 
to the alkyl group. The large electronegativity of fluor 
causes a reversal of the dipole in CF3S and CF3CH2S. 

We deflne the contribution to the interface dipole re- 
sulting from chemisorption as /ichcm ~ A/i — /^sam ■ The 
results shown in Table ITTl clearly demonstrate that /ichcm 
is nearly independent of the molecule and strongly de- 
pendent on the metal substrate. As an independent 
check we have also calculated the dipole on the basis of 
the electron density redistribution, see Fig. ^ fj,^„ = 

~^ IIIcoii^^'^(^)'^^'^y^^- We flnd that fi^n ~ A^chcm, 
which indicates the consistency of this analysis. 

The results obtained allow for a simple qualitative pic- 
ture. The chemisorption dipole /ichcm is very small for 
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FIG. 2: (Color online) Work functions of clean surfaces Wcican 
(blue), including the chemisorption dipole Wchcm (black), and 
of the SAM covered surfaces (red, green). 



all SAMs on Au(lll), indicating that the charge transfer 
between the Au surfaces and the molecules is small. This 
generalizes previous results obtained for methyl thiolate 
SAMs on Au(lll) [11 Q. Since the work function of 
Ag(lll) is substantially lower than that of Au(lll), a 
significant electron transfer takes place from the surface 
to the molecules for SAMs on Ag. This is confirmed by 
the values of ^chcm for Ag in Table ITTl Fig. ^ shows that 
the electrons are transferred mainly to the sulphur atoms. 
Integrating the positive peak of An on the sulphur atom 
gives a charge of (— 0.24± 0.02)e. The sign of the charge 
transfer is such that /.tchcm increases the work function 
with respect to clean Ag(lll). By a similar argument, 
since the work function of Pt(lll) is much higher than 
that of Au(lll), an electron transfer takes place from 
the molecules to the surface for adsorption on Pt. The 
values of /ichom for Pt in Table Hll confirm this. In this 
case the net charge on the sulphur atom is positive and 
A'chcm decreases the work function with respect to clean 
Pt(lll). 

The size of the charge transfer is remarkable. 
Chemisorption creates an interface dipole /ichom that 
overcompensates for the difference between the metal 
work functions. We define a work function that in- 
cludes the contribution from the chemisorption dipoles 
as Wchem = Wciean + e/Zchem/ (eq^) ■ The rcsults showu 
in Fig. 13 demonstrate that Wchem decreases in the order 
Ag, Au and Pt, whereas Wcioan increases in that order. 
The work function of the SAM covered surfaces can then 
be expressed as W = W^chom + e/isAM/(eo^)- From the 
polarity of the molecules discussed above, it is clear that 
SAMs of CHgS and CH3CH2S decrease the work func- 
tion, whereas SAMs of CF3S and CF3CH2S increase it. 
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tions and B. de Boer, P. W. M. Blom and P. J. Kelly for 
very helpful discussions. This work is part of the research 
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